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B | Abstract: The third design procedure for fully cavitating hydrofoils 
is based upon a linearized inverse theory of two-dimensional 
cavity flows at arbitrary cavitation number. The cavity 
surfaces are assumed to originate at the leading and 
trailing edges of the wetted surface. This report reviews 
and completes the basic theory and gives detailed examples 
obtained from the resulting parametric design technique. In 
the third design procedure, one specifies the design lift : 
coefficient, the cavitation number and the upper cavity 
thickness at two points along the profile chord. A prescribed 
pressure distribution shape is also selected. These 
quantities determine the profile design which consists of 
the upper cavity and wetted surface contours, the design 
angle of attack, the cavity length, the drag coefficient, 
the moment coefficient and the lift-to-drag ratio. The 

{ method also includes off-design calculations in accordance 
with the direct linearized theory of cavity flows which determines 
the flow states for which interference can occur between the 
upper surface of the cavity and the upper nonwetted surface 
of the profile. The hydrodynamic performance of specific 
“point designs" is also given by these direct calculations. 


The chief new feature of the third design procedure is the 
ability of the designer to prescribe two points on the cavity 
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surface instead of one as heretofore. Although certain 
constraints must be observed by the designer when specifying 
these two values of cavity thickness, the third procedure is 
found to be more general and more flexible than the first or 
second procedures studied previously. The necessary constraints 
are incorporated in the computer logic for the method. The 
fact that linearized theory is used tends to limit the 
applicability of the procedure to conceptual design and 
feasability studies. The computer program for the procedure 
has been found to be economical and well suited for its 
intended purpose. 


This research was carried out under the Naval Sea Systems 
Command General Hydromechanics Research Program, Subproject 

SR 023 01 01, administered by the David W. Taylor Naval 

Ship Research and Development Center, Contract N00017-73-C-1418. 
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Nomenclature 


cavity length parameter 


section lift coefficient 

section drag coefficient 

section moment coefficient 
cavitation number = (P|-P )/(pu2/2) 
cavity length, & = ie +1 


lift-to-drag ratio 


pressure amplification factor used in the second design 
method 


nondimensional perturbation pressure 

cavity pressure 

free-stream static pressure 

peak pressure location, measured from the nose 
upper (non-wetted) surface of the hydrofoil 
cavity thickness at the trailing edge 

free stream velocity 


function used to compute the angle of attack corresponding 
to the cavity-foil interference in the off-design calculations 


absissa, taken as positive downstream, with the origin at 
the nose 


distance of the center of pressure from the nose 
chord wise location of second cavity-thickness tontrol point 
cavity contour (upper surface) 


cavity thickness due to the complementary function (not 
total thickness) 


complex variable in physical plane, z=xtiy 
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Greek Symbols 


angle of attack of the chord line with respect to the 
free stream 


cavity parameter, l - =. 


V1+K 
wetted surface shape 
ratio of cavity thickness to T at X=X 


density 


3) (a,x) ; 


5, (a,x) functions used to compute the cavity contour in the off- 
design calculations 


Superscript 


. prime indicates differentiation with respect to x 
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INTRODUCTION 


This study is the second of two investigations intended to provide 
a two-dimensional linearized theory and companion numerical methods [1] 
for the preliminary design of supercavitating or superventilated hydrofoil 
sections which will be consistan: with the linearized mixed-foil theory 
of Wang and Shen [2,3]. Reference 2 describes the design of hydrofoil 


sections which embody good performance in both fully cavitating and 


noncavitating flows. In order to provide a profile of high lift-to-drag 


ratio in a cavity flow the designer prescribes the shape of the lower 
surface with the foreknowledge that his prescription of the wetted 
surface contour will achieve the intended performance. On the other 
hand, he must provide a low-drag section in noncavitating flow and also 
provide a profile which can tolerate a prescribed sea-state without 
suffering cavitation inception until after the take-off speed has been 
reached. Thus the pressure distribution of the upper surface of the 
noncavitating hydrofoil must also be controlled by the designer. 

The linearized mixed-foil theory of Wang and Shen [2] pertains to 
two-dimensional foils in an unbounded fluid. The lower surface contour 
is specified in terms of high-speed superventilating or supercavitating 
performance and the upper surface pressure distribution is specified in 
terms of sea-state requirements for moderate speeds. The hydrofoil section 
of streamlined shape is then computed from the theory. The purpose of 
the present investigation is to arm the designer with the prerequisite 
knowledge of wetted surface shape and related cavity-flow parameters 
so that he will have a chance to execute his preliminary mixed-foil design 


without inordinant travail. 


: 


= 
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The first of the two studies noted above was a numerical design study 
[4] which explored the hydrodynamic aspects of fully cavitating hydrofoil 
section design in accordance with the linearized inverse theory of cavity 
flows [5,6]. In these calculations the cavity springs from the leading 
and the trailing edges of the hydrofoil so that the upper surface of the 
foil is contained inside the upper cavity contour and only the lower surface 
of the profile is wetted by the flow. A typical flow geometry and a 
prescribed pressure distribution is illustrated schematically in Figure l. 
Two design procedures were used in Reference 4. In the first of 
these procedures, the design is carried out at the ideal attack angle 
(or at "shockless entry"). The thin-airfoil analog of this first design 
procedure is given in Reference 7, In the cavity-flow problem the airfoil 
thickness distribution is replaced by the cavity thickness. In the first 
design procedure interference between the upper surface of the hydrofoil 
and the upper cavity contour is eliminated by use of the "point-drag" 
solution of linearized cavity-flow theory. In the second design procedure 
cavity interference is eliminated by increasing the attack angle above 
the shockless entry condition. The point-drag solution is not used in this 
procedure. In both of these design procedures one prescribes the cavity 
thickness, T, at the trailing edge, the cavitation number, K, the design 
lift coefficient, Ch» and the pressure distribution on the wetted surface 
of the profile. Then the hydrodynamic moment coefficient, Cu the cavity 
drag coefficient Ch and the section lift-to-drag ratio are calculated. 
The design attack angle, the wetted surface contour, the upper cavity 
contour and the cavity length also follow from the theory. 
The theoretical developments of References 5 and 6 also contain the 


ingredients of a third design procedure in which the cavity thickness is 
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controlled by a combination of the point-drag solution and the use 
of an attack angle which is greater than the ideal angle. Evidently 
the third design procedure contains the first and second procedures as 
special cases. Until now the third procedure has not been worked out 
in detail even though, as we shall see below, this procedure offers the 
designer somewhat more control of the upper cavity contour than he can 
obtain with the first or the second design procedures. Therefore, it 
is useful to complete the theoretical development for this procedure and 
to program the resulting numerical methods for a computer [1]. In the 
following sections, this development is described. Then some hydro- 
dynamic and geometric consequences of the third design procedure are 


presented. 


THEORETICAL PRELIMINARIES 

In order to give a reasonably self-contained account of the third 
design procedure, we will outline the basic theory which underlies the 
method. Further details are given in References 5 and 6. 
Flow Geometry and Characteristic Parameters 

In the present linearized theory, we work with velocity perturbations 
created by the Palbeonwier system in an otherwise steady rectilinear 
flow. Let us suppose that the chord of the profile is inclined at the 


small angle a with respect to the free stream velocity, U, as illustrated 


in Figure 1. The flow direction is from left to right as shown. Let 
us fix an x-y coordinate system at the profile nose and with the x axis 
parallel to the free stream direction. In this coordinate system, the 


trailing edge of the wetted surface (of unit chord) will be at x=1l, y=-a. 
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With respect to this system, we can write the free stream velocity 


as 


> 
q,, = {U,o} 


The magnitude of the velocity on the cavity surface is 


ie UV1+K 


which follows from the definition of the cavitation number, 


and the application of Bernoulli's equation with P,. being the free stream 
Loe} 


static pressure, P, the cavity pressure and p the liquid density. 


k 
If we use the cavity speed q. to normalize the velocities we can 


write 


1+u 


a(x,y) ics 
nv, 


where the dimensionless disturbance velocities u and v are created by 


the foil and its cavity. Far upstream of the foil, the perturbation velocity 


components are 
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ie 


V1+K 


Let us define the dimensionless perturbation pressure p as 


: P - Pe 
le 


1 
zea 


Then to first order terms, we must have 


p(x,y) = -2(1+K) u(x,y) 


and on the cavity, p=u=0. 
Conformal Transformations and the Complete Solution 

We can use certain conformal mappings in order to transform the domain 
of the flow in the z-plane (z=xtiy), in such a way that many key elements 


of the solution for the complex disturbance velocity, 


w(z) = u-iv (7) 


can be found by inspection. In these transformations the complex velocity, 
w(z), is taken to be invariant at corresponding points. The method of 
solution is conveniently illustrated by a study of the geometry of the 
flow boundaries in the various planes. The first step in this process is 
the representation of the flow boundaries in the z-plane in a way that is 
consistant with the linearization implied by the use of the perturbation 


quantities u, v and p. All of these will be much smaller than unity as 


-16- 5 July 1977 
BRP: JF:jep 
long as a and K are also small quantities. Generally the wetted surface 
of the profile will be cambered with a camber distribution n(x). The 
magnitude of n will also be much less than unity at every x along the 
wetted surface and in particular, n=0 at x=0 and x=l. That is, n is 
measured with respect to the profile chord line. 
In order to simplify the following considerations let us imagine 
that the total complex perturbation velocity consists of two parts: the 
first part being that of a flat plate profile at an attack angle ao and 
the second part being due only to the camber distribution n(x). To first 
order accuracy we can replace the foil and cavity in the x-y plane of 
Figure 1 by a cut of length & along the real axis as shown in Figure 2. 
This figure also shows the boundary values which apply to w(z) for the 
flat-plate part of the solution. We can transform the flow outside the 
cut in the z-plane into the region outside of the unit semicircle in the 
upper half of the f plane. As illustrated in Figure 3, the transformation 
can be carried out by a sequence of mappings. In the first of these, we 
transform the point x=% into the point at + with a bilinear transformation 
which leaves the points 0 and 1 invariant. Then we take the square root 
of this mapping to get the configuration shown in the V plane of Figure 3. 


One can verify that the resultant mapping is 


Since the foil and cavity are mapped into the real axis in the y plane, one 


now applies the Joukowski transformation, 
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in order to map the wetted surface of the foil into the upper unit 
semicircle in the f plane. The boundary values on w which would apply 
to a flat plate hydrofoil at corresponding points in the z, V and [ 
planes are so rte in Figures 2 and 3. 


In the f plane, one can write the solution for w(t) as follows: 


me ae See! Ae a n 
w%) =7y tig & z) + i(ata) i aera ese} 


In order to see that this is the solution for a cambered supercavitating 


profile, we note that the first term, 2iA/(T-1), is equal to 


A[cot o i] 


on ewe’, Therefore its imaginary part is simply the constant -A on the 
unit circle. On the real axis where [ is real this term is purely imaginary. 
Therefore it satisfies the condition u=0 on the cavity. “The second term, 


iB(t-1/0)/4, equals 


on the unit circle. Therefore it does not contribute to the imaginary 
part of w(t) on the wetted surface of the profile. Moreover, when % 

is real this term is purely imaginary and so it also satisfies the 
condition u=0 on the cavity, as does the third term and also the 

Laurent series, provided that the C, are real. The first three terms from 
Equation (10), and the Fourier coefficient Cy can be thought of as 


those terms which allow for the angle of attack (or flat plate) part 
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of the solution. We note also that these terms also provide for the 
branching of the streamlines at the profile nose C=1 and at the end of 
the cavity, x=£+>T=~, Both of these points are singularities of w(Z), 
as one expects for linearized stagnation points. Thus the first term 
involving the constant A opens up the cavity and the second term involving 
the constant B closes it again. The Fourier series represents the contri- 


bution of camber to the solution. In the case of the direct problem 


the camber function n(x) is known in advance. The coefficients C. are 


determined by the derivative dn/dx in the usual way. 

Of course, in the present inverse problem the camber function n(x) 
remains to be determined. Instead of specifying n, we prescribe the pressure 
distribution p(x) on the foil. Therefore we will modify the form of 
Eq. (10) to make the specification convenient. This modification is guided 
by the fact that on the wetted surface of the foil the expression (12) shows 
that the second term of Equation (10) contributes to u but not to v when 
eee, Therefore, if we add a term iB/2t to the solution, the above 


contribution to u(@) will be cancelled by the added term and we can relate 


the prescribed pressure distribution, 


s 
, 


p(9) = -2(1+K) u(6) 


to the camber term in Equation (10). Note that we will still have a 
separate contribution to the pressure on the profile from the first 
term of Equation (10) which is due to attack angle as shown by the real 
part of expression (11). This part of the total pressure distribution 


can not quite be prescribed arbitrarily although one can find an ideal 
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attack angle at which this first term (and the resulting singularity 
at the profile nose) will vanish, corresponding to shockless entry. 
Finally we do not wish to be restricted to a Fourier series represen- 
tation for the pressure distribution p(9). Therefore we will write the 
modified complex velocity as 
2iA 


w, (©) = ia + O48 - ee + iD + w, (2) 


where 


The term w, (2) is an analytic function of f which is regular everywhere 

in the flow and represents the effect of profile shape or pressure 
distribution upon the particular solution 5" Now the particular solution 
has been designed to satisfy the boundary conditions on the profile and 

on the cavity as discussed above. We can add to this particular solution 
any other solution which is regular at infinity, which satisfies u=0 on 
the cavity and on the foil, and which produces a closed body in the 


interval 0 <¢ x < &. Such a complementary function is given by 
B ptf 2 2-2 
a mer PP 
wz) = - 5 OVSy ys (14) 


and it is known as the point-drag solution of linearized cavity flow 


theory. The complete solution of the problem is then given by the sum 


of Equations (13) and (14): 
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w(z) = M, re Wo (15) 


Boundary Conditions for the Inverse Problem 


Equation (15) is the form of solution which must now be applied so 


as to satisfy the boundary conditions on the present inverse problem. 


r These conditions can be listed as follows: 

. (i) w(z) = - 8 at z= 

L (ii) u(x,0+)=0 on the cavity surfaces | 

" (iii) v=-at a on the wetted surface of the foil ‘ 

E (iv) the cavity together with the wetted surface of | 
; 


the foil is a closed body. Thus 


In § w(z)dz = ~2Ma, = 0 ; 


where ay is the real part of the residue of 
i w(z) at z=0. ‘ ' 
(v) w(z) is continuous at z=1 (Kutta condition). 
f (vi) us<0O everywhere in the flow. In particular, the 
cavity pressure is the lowest pressure in the 
flow. § 
| Before we can apply these conditions and determine the various constants 
such as A, B, D, and E in Equation (15), we must determine the analytic 


function w, (2) in terms of the prescribed pressure distribution p(x). For 


our purposes, it is sufficient to determine the value of w, (2) on the 


cavity and on the foil. In the y plane, this can be done by means of 


the Hilbert transform or by certain direct integral-superposition 


a lial 


=e 
” 


es og 


¥ Sait, vars aaa 


ag ep 
= P 


— ES 


in the integral superposition without calculating w 
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techniques as were used in References 5 and 6. In either case, the above 

boundary conditions give a system of equations involving definite integrals 

in the vy plane. These equations can then be transformed back to the 

z plane with the hale of Equation (8). Because we are working with the 

function Wy only, we can work with the various boundary conditions directly 

1 explicity. However, 

we can not apply integral superposition to the problem as a whole because 

of the nonlinear involvement of the cavity length in the problem. 

In any event we obtain from conditions (i) chrouph Gri) a series of 
equations involving integrals over the wetted surface which contain the 
function p(x) as one factor in the integrands. We will simply list the 
results below; but before we can do so we need to determine the relative 
importance of the individual contributions from the prescribed pressure 
distribution and from the angle-of-attack or flat-plate part of the 
solution to the design lift coefficient C,. In order to evaluate these 
effects, let us denote the contribution of the prescribed pressure 


distribution by 


1 
Cr = | p (x)dx , (16) 
0 


Then the difference, Cr - Crs is equal to the contribution of the flat 
plate part of the solution and it is known that this contribution will 
be equal to 4m (1+K)b, » where by is equal to the imaginary part of the 


residue at z=0 as derived from Equation (13) with the term Wy deleted. 


The point-drag solution does not contribute to the lift. ‘The final 


result can be written as 
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(i-mc, = - 4 n(x) S23, (17) | 
a 
where 
mC, = ct (18) 
and 
9 S . aiae 1/2 
= [Y2 + 1] + [V2 - 1] ‘ (19) 
€ 


Returning to the employment of the boundary conditions, one finds that 


condition (iv) leads to 


a 
pe - fete, a | Cr a, 


Q . (20) 
4 alva 47 (1+K) a& } Jets) 3 


provided that p(0)=0. After some manipulation the two equations resulting 


from condition (i) can be written as 


L 
B fee et Briain Bee 
ge ie 2 Pane | > vet? 
0 
and 
(1-m)¢; (2a743) tad m f R-x 
(1+K)86 = oo jane + (14+K)E + ont | ry p(x) dx r (22) 
0 


Condition (ii) has been satisfied by the construction of the solution, 


Equation (15). Moreover, condition (v) is satisfied by the flat plate 


part of the solution as can be seen from Equations (11) and (12). It will 
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also be satisfied by that part due to the prescribed pressure as long 
as we require that p(1)=0. In fact it has been found [5] that the strict 


condition to be satisfied is 


% péx) ~ 204") tite (23) 


near x=l, where A>0 may be selected arbitrarily. Condition (iii) gives the 


equation for the shape of the wetted surface. The result is 


Ps | 
+ — 1. 
: a oe ha fei ci ny a ee) 
n(x) = (C, - 5)x + aB[% tan Vee ¥x(2-x)] - Ev¥x(2-x) + 2n(1tK) }2 | Cox 
2 a‘! 1 
q -} (7x (2-x) + Qtan * Nf el | P(x) dx + ane \ = | p(x) Oi eee dx 
9 Vx (=x) 2 ‘ 
c 1 = 
_L Ey sy /x(&-t) 
+3 in) | p(t) &n |1 “ee dt , (24) 
0 | 
where we have used Equations (16) and (18) in order to introduce the term 
which depends on Ce The satisfaction of condition (vi) will be strict 
in the following design method. We shall insist that the contributions 


of both flat-plate and prescribed pressure distributions shall always exceed 

zero on the wetted surface between the leading and trailing edges, 0<x<l. 
One last formula which is an important ingredient in profile design 

is the upper cavity contour. It is obtained from considerations similar 


to those leading to n{x). When the cavity ordinate is measured from the 


profile chord line the result is 
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y(x) = ax - 2a” RAF, (x,a) + aBl[Yx(£-x) - re ! x Ca + D)x 


1 5 ace. ; 
+ E yx(2-x) + ETERS) G(x,a;p) 


1 
G(x,a;p) = a | BG) ds +5 [vx(2-x) + iia > 


-| pct) Sa 11 4 ees dt 
0 


saay 8 B 
Fy(x,a) = yp POD + —F Levey + 4 10,29] 
2a V'ta 4w 2w 


with v=a Ve for x20 on the upper cavity surface and w? = ave, 


By =~) / daa} /2 ; By =- -\/ (w24+a")/2. The L and T functions are 


2 {y,v + 6, - ACH +¢y,v + 6)” 
L(v,a) = &£n 4a [——————- 5 REC 5} ] 
(55 + 6) +a) 


T(V,a) = tan - 
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oe ab, + B,/2 et so aB, - B,/2 
‘ae 2 9 


THE DESIGN PROCEDURE 

As explained at the outset, the third design procedure contains the 
first and second procedures as special cases because both flat plate and 
point drag solutions are used to control the cavity geometry. Thus if it 
should happen that A=0 or E=0, these special cases will be the outcome 
of calculations based upon the prescribed input parameters. Neither A nor 
E will have been set equal to zero from the start. The various parameters 
which have been identified in the problem thus far are the constants A, B, 


Cy» E, a, 2, K, m, C. and the functions p(x), n(x) and y(x). 


L 


In the formulation of the design process the two constants K and Cc. 
will be prescribed from the outset and of the three functions, p, n and y, 
only p(x) will be an input quantity. The outputs of the design process 


will be the geometric parameters, a and & and the cavity and profile shapes, 


y(x) and n(x). The hydrodynamic performance parameters Ch» ce and related 


parameters such as L/D and center of pressure location x at the design point 
will also be considered as primary output quantities. 
Provisions for a Determinate Procedure 

Naturally, the determination of the above primary output quantities 
depends upon the evaluation of A, B, Cy, m and E in terms of the 
inputs. Thus, when we add these five parameters to the two unknowns, a 
and m, we find that there is a total of seven unknown quantities to be 
determined. At present, we have found only four relationships between 


them. These are Equations (17), (20), (21) and (22). Three more relation- 
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ships are needed to produce a determinate design procedure. These may 
be found from additional conditions on the functions n(x) and y(x). 

For example, consider the camber function n(x). We have stated 
that these ordinates are to be measured from the profile chord line. 
This can only be true if the arbitrary constants in Eq. (23) satisfy 
explicitly the condition n(1)=0. Thus, we must have the added relation- 


ship, 


O=a- (G +d) + able tan? = - a] ah + Sey 


1 1 
C 
- F(a +2 cau 4) | —P(x)dx 5 tan 2 | p(x) 7. + = tn 4 
0 vx (L-x) 0 a 


1 
. | ote) de fl + 4 /=* a (30) 
0 


Another point is that we have not yet made any provision to insure that 
the upper cavity contour will clear the wetted surface of the profile 
and still leave enough room for the hydrofoil structure inside the cavity. 
Since we require two more relationships to make the problem determinate, 
we can impose two conditions on the upper contour of the cavity and use 
them to control the cavity clearance. One of these can be the specification 
of the cavity thickness at the trailing edge of the profile. Recalling 
from Eq. (25) that y(x) is measured from the profile chord, we take y(1)=T, 


where the cavity thickness T is now an input parameter as it was in the 


first and second design procedures [4]. The resulting relationship is 


35 tm te 
; 
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(1-m)C 2 (1-m)C 
fu L 4a” va - 1 L e+ad 
Peas m(1+K) ° 6-ae F)G,a) + [ark tan a )2m(1+K) ° £(6-ae) 
1 (1-m)C 1 
ue | Bete) | 8} ea 
4m (1+K) 2 Je w(1+K) * 6-ae 4m (1+K) L-x 
3 0 0 
(1-m)¢, a (2a743)¢ + ad m a ; Qa 
Wee Gees 8 eee ~ Or(tk) * al Ug Biase 
0 
m 
+ SrqaKy CG.asP) (31) 


where the functions Fy and G have already been defined in Equations (26) 
and (27). 

Next we will require the cavity to have a specified ordinate equal 
to wT, measured from the profile chord line at a point Xo Near the nose. 
Obviously the factor vw will be in the range O<p<1l. We can not specify 
the value of with complete arbitrariness. For example, if the cavity 
pressure is to be the lowest pressure in the flow field, it follows that 
the cavity must be convex. Therefore if T is prescribed, it follows 
that we must have U?X 5» because HEX) corresponds to y(x,) being on a 
line drawn between the nose and the point (1,T), where these ordinates 
are measured with respect to the chord line. As a practical matter, it 
is desirable to make the cavity as thick as possible in the vicinity of 
the nose so that the leading edge of the profile can be strengthened. 


The condition y (x, )=HT is certainly not the only possibility for the 


last condition needed to make the design procedure determinate. For example 


(1-m)C 2, 7 F,(x_,a) (2-x_) ' 
Ut. L 4a°&/a “1°"0 ae se. | 
anil T(1+K) ° 6-ae Xo is x x, “eh \ 
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an alternative condition which we have not investigated but which could have 
some advantages in a design procedure, is y (x )-n(x,)=uT, where w is a new 
parameter which measures the clearance between the upper surface of the cavity 
and the wetted surface of the foil at the point Xo" This alternate condition 
could be used to make certain from the outset that there is no interference 
between cavity and foil near the nose. The condition which we have adopted 
in the present method will not permit us tc guarantee non-interference from 
the start. However, as we shall see, UW can be prescribed in a limited 
range which depends in each instance upon the input parameters K, Cc and 
T as well as upon the shape of the pressure distribution, p(x). Therefore, 
one éan avoid interference by designing a number of profiles for various 
U-values in the permissible range and selecting the preferred design from 
the several possibilities. 


Use of the condition y (x,)=uT in Equation (25) leads to 


oO oO 
(1-m)C , (1-m)C 
saci? e+ad + m (2-2x) p (x) q . mre 
Qn(1+K) * &(6-ac)  4n(4K)R Eat eT HIaR) | O=8e 
0 
1 = (1-m)C 2 
a st. | Biedde fo Jy . Ga tapered _ 
4m (1+K) Q-x Xo 27 (1+K) 2 6-ae 
0 
1 
G(x_,a;p) 
pee Oey ee Ritts sei 
~ 2m (14K)2 | < Peery * x, 2m(14K) ; ae 


In Equations (31) and (32) the parameter E has been eliminated by using 


Eq. (22). The same step could also have been taken in order to eliminate 


len i ote A 
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E from Equation (30). Equations (30), (31) and (32), when combined with 
Equations (17), (20), (21) and (22), provide a determinate system. 
The Sequence of Calculations 


The preceding formulae permit the determination of all unknowns for 


F a profile design. This determination can be shown to depend solely on 

the cavity length parameter, a, by means of a tedious process which 

a eliminates the other unknowns and results in a nonlinear equation defining 
3 a- In an important step in deriving this key result, one eliminates a 


from the system of equations by subtracting Equation (32) from Equation 

? (31). Then Equations (17), (20), (21) and (22) can be used to reduce this 
equation to one involving only the unknowns m and a. Operating with these 

¢ first four equations on Equation (30), one can reduce it to a second 
equation involving m and a Fortunately, m can be separated out as a 
linear factor in the combined equations resulting from the difference of 
(32) and (31). The reduced form of Equation (30) can be factored in the 
same way. One then eliminates the quantity m by dividing one of these 
equations by the other in order to obtain a single nonlinear equation with 
the only unknown being the cavity-length parameter a. The resulting 


expression is: 


Q-x Cc 2 F_ (x _,a) 
-~ Hyp _ af Oi enti. ee Rva ee pet 
me! — (a Xe )B tT (1+K) 6-ae (F, (a) ) 


x 
co) 

+ [a - 2 ao = Q(tan t Ew b. tan? \/ “o )] Stns. 
x aX, dax 22(6-ae) 


5 July 1977 

BRP: JF: jep 
Lx, C. - 2 
ie (2a (2a743)e+ad fi. «ft < Sees L 4a“ Va 


we “Sitisaey T(+K) }~ écae Fy (+a) 


+ 


2 
a (2a +3)etad hb Ly etas 
aaa farhten” 2 Tiscat! Be 


3 
2 -1 1. e+aéd a (2a +3)et+ad 
F, (1,a) + (a - 2 tan ) R(d-ae) ~ d-ae 


(2-2x) p (x) an sé 
Vx (2-x) 


1 
-F (a + & tan +) [ Boose ete 


0 ¥x (2-x) 


avx + /R-x 
avx - VR-x 


R-2x 


p(x)dx + = 
0 ¥x(2=x) 
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and where the functions F, and G are defined by Equations (26) and (27). 
The known quantities in Equation (33) are K, Ch» Te ts X and p(x). The 
unknown quantity is ing a0. The value of a can be found by trial and 
error for each set of known quantities. 
We can also introduce the additional notation 
o ha’ tva et+aéd 


- _“L_ |4a‘tVa 2 mh Boy «Se 
Hy tee eee) feeacy Teed + = * Bee”? gee) 


(2a°+3)e+ad 


_2 
g (6-ae) 


in order to write the simultaneous conditions n(1)=0 and y(1)=T from 


Equations (30) and (31) as 


a -m 
Fy ask) My 


Similarly, we can introduce the expression 
R-x C 2 
“ a Bete 2 as? oy. L ha iva 
~~ Si Bla x. ) ~ 5G) | bac 


etad 
22 (6-ae) 


= (2a743) e+a5 (a - 
22(6-ae) 


(35a) 


and write the simultaneous conditions y (x) = uT and y(1) = T from Equations 


(31) and (32) as 
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—m 
Fy 1 (1+K) Hy ; 


Then in accordance with the discussion at the start of this section, we 
can divide these two expressions in order to eliminate m and write 


Equation (33) in abreviated form as 


f(a) = F,H - FH, = 0 . (33a) 


Once the solution of Equation (33) has been found, we can solve for 


m from the condition following Equation (35): 


a eo (36) 


Then in accordance with Equation (17), we can find A from 


~ “W(4K)  §-ae 


The design attack angle follows from 


(1-m)C 2 
a L 6 : a -1 1, etad _ a_ (2a +3)etad 
me ae eee) Yeon aa  * Pet 9D) Sege OE 
1 1 
m a IQax ui =I R-2x : 
+ In(1ak) )7 2 | VV = p(x)dx + a) (a - & tan _ | Mise p(x)dx 
0 0 
c 1 1 
a alt Qn _ Lis + 2 tan’ dy —P(x)dx _ tan” = p(x) —~ dx 
2 2 a a R=x 
a 0 x (L-x) 0 


1 
+ | p(x) &n (1 - Ivo ‘ (37) 
0 
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After we have calculated m, a and A, we can calculate the remaining 


parameters in accordance with 


1 
A etad m | (2-2x) 
een a ees See “ p(x)dx , (38) 
4 an fe 4m (1+K) ak Jeticay 
1 
ents te cir ee ee ae a Bee 
ee. te 2 + Tra | i ox (39) 
* 0 
and for E = y,(1)/a we have 
semiecbe (2a7+3)e+a8 m i Sox 
i= 8- aay ey ee | = a 
0 


The hydrodynamic force coefficients C, and Ch follow from 


D 


(1-m)C 
Cy = 27(1+K)2 (aB + 5? = aay + cia + ae | wo oe p(x)dx 
2 Vx (R-x) | 
| 
+ T(14K)E ’ (42) 
and 
2 2 ) 
c. = ~(1-m)¢, (3ta")6 - 2a(2ta je - mm | xp (x) dx . (42) 
8ev2 0 


The wetted surface contour, n(x), and the shape of the upper surface of 


the cavity, y(x), follow from Equations (24) and (25) respectively. In 


addition, the slope of the wetted surface is 


> sine | —- = ve 


wr} 
> 


oO 


1 hatte ode 


1 
ee = Eo = & see m 1 | p(t)dt 
Bg ee GD <8 ot ae | tt 


-34- 5 July 1977 
BRP: JF: jep 


vx (i-x) 0 \q 
1 , C. x 1 
- $= | ee 2 f pce) \gte at + ay i 
* f veti-t) = ave(imx) f = neem? 
1 
+ Beek) ee Pp | pee) ae ‘ (43) 


avx (=x) 3/2 , “\e- ~ zs 
=t =x 


where P denotes the Cauchy principal value of the integral so designated. 
This result is required for off-design calculations of profile performance. 
For reference, one can also tabulate the perturbation pressure on the wetted 
surface. This pressure is the sum of the prescribed pressure m p(x) and 
the pressure p(x) which is due to the nose singularity. This latter 
pressure is obtained from those terms involving A in Equation (13) when 


tae? or when wine, The expression for P(x) is 


aehisRlides R-x ~ avx K>0 
1 (6-ac) di MY 
p(x) = 
4(1-m) C 1-Vz 
ui vx ’ =0 ’ 


where the second form for p(x) is simply the limit of the first as K*0 and 


as a’, Then the total pressure is p(x) 4+ mp(x). 


alia te ee 


—— 
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Formulae for Profile Design at Zero Cavitation Number 


The results appropriate for this special case can be obtained directly 


from the preceding formulae by taking the limits K>0 and 2*~. The 
simplifications introduced by this step indicate that when K=0, all 
numerical work should be based upon these specific formulae. 

For example, it can be seen from Equation (20) that B?O as K>0, and 
since a’ in this case no iterative solution such as that associated with 
(33) or (33a) is needed. Therefore we can proceed at once to the 
calculation of m: 
oe 2c, ae - ua 2u/x IT 


2 C.J (x )+ oe - 1) | p(x) &n bas dx + | p (x) [inti - i &n(1+/x /x) ]dx 
fe) 2 = Je x fe) 
fe) 0 


1-/x 0 


(44) 
where 
J(x,) = ie! + Ly [3/2 = en(1 + 72)} + “a an (xen + Vit) 
x (e) 
oO 


«~ 4 2vx 5) Vi + Re [3 . (45) 


Then we can find the design attack angle from 


1 1 
a= y (1) + 2Gem o 4 B- jf ROD 4 | ox) fn t~Segeh (46) 
a ai 0 = 


We have noted above that a*™ as K*0. This fact requires that we return to 


the meaning of the coefficient E in terms of its contribution to the cavity 


thickness, T=y(1). As we have seen in Equation (25), the contribution of 
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the point drag singularity is Evx (2-x) = aE at x=l. If we denote this 
contribution to T by the designation y,Q@), we can write E=y (1)/a. That i 
is, recognizing that y,@2) will be bounded as K*0, we require that E=0(1/a) | 
in this limit. We note that this finding is consistant with the result 
which can be deduced from Equation (22), for example. Now if we replace 


E by y,Q)/a in Equation (25) and let a*~ and x=1 so that y(1)=T we get 


al 1 
(1-m)Cc : = 
yh) = 5 - a [372 - na + 72)) - [ RGdd 2] pis) tn SEE aah . cans 
0 Vx 2: 0 1-vx 
Also in this case, we see that 
2(1-m)C, 
A=- 7 ° (48) 
Then we can write 
2(l-m)c, : aids ; aa m 
n(x) = (a - ee - ao ive | a. + | p(t) &n |1 - vx/t}dt} - yD) ¥x (49) 
a 0 


and 


(1-m)C cr 
y(x) = ax + —— {(1427x)x2/4 Vi + va = 2x ~ teat? + Vi + vx )} 


1 1 
+ y (1) vx + ~ (| P(x)dx _ | p(t) &n (1 + Yx/t)dt} . (50) 
0 vx 0 


The hydrodynamic coefficients C. and Ch. become 


D 


D 


1 
C. = + {2(1-m)C, + a | AR dx + ny (1)? (51) 
0 


vx 


= 
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and 
1 
SS oe hie: 
C. = 16 (1 m)C, m | x p(x)dx ‘ (52) 
0 
Note also that the slope of the wetted surface is 
2(1~m) 9) ng f (tat f (x)dx 
EC Rm Sa ee | se a | saps ae (53) 
2vx 4nvx 0 vt - vx A vx 


Preliminary Analyses 


The analyses of the preceding two sections have been formulated with 
the thought that one specifies the design pressure distribution, p(x), the 
cavitation number, K, the cavity thickness at the trailing edge, T, and 
the design lift coefficient, Ch. He also prescribes the cavity thickness, 
uT, at the chord location, Xo" The computation sequence outlined above then 
produces the cavity length, &, the wetted surface shape of the hydrofoil, 
n(x), the upper contour of the cavity, y(x), the design attack angle, a, 
and the hydrodynamic performance parameters Cy Cy and L/D. 

However, one is not able to specify the parameter uw with complete 
arbitrariness. This situation arises as long as one insists that the 
pressure on the wetted surface is greater than the cavity pressure except 
at the trailing edge, and possibly at the nose. The wetted surface 
pressure will certainly be positive if the parameters y.Q@) and m satisfy 


the conditions 


y,(1)20 ’ 


(54) 
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Depending upon the specific values of C_, T, K, X and the nature of p(x), 


tr? 
the added constraints specified by (54) indicate that permissible 


values of uw will lie in restricted ranges. Therefore as a prelude to 
the execution of any profile design, the designer needs to have some 
knowledge of the permissible range of ) values which he can specify. 
The purpose of this paragraph is to explore this matter in detail 
sufficient for us to continue our study of the third design procedure. 

-The basic nature of the constraints noted above is probably most 
conveniently illustrated when K=0. Then we can determine the relationships 
between Cy» T, y,@), m, X, and uw quite simply. In fact the basic equations 
needed to carry out this investigation are given by Equations (47) and (50), 
with the latter of these being written for the particular case in which 
y(x,)=HT. Equation (46) is used to eliminate a from the final expression. 


As a result we find that 


¢ — 
1 fe hie 372 - &n(1t+v/2) _ _m_ 


ie 1 


Cc plccaaiaiine 
+ oe JAB (arte) Vx otva, - x, [372-2n (1472) ] 


1 1 
- taxi! 4 Vi1tvx°)} + { (Vx -x,) | — - | p(x) 
L x 
0 0 
<i a 
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C 
Thus when K=0, we can prescribe values of m, = and X and determine 


corresponding values of y,(1)/T and Uw. 
bs Rectangular Pressure Distributions 
In order to continue the analysis we can introduce a specific pressure 


distribution into Eqs. (55) and (56). For our purposes, it is particularly 


convenient to study pressure distributions of rectangular shape. We will 


in specify these by the following equations: 


(a) for O<s<1/2, 
t C,/2s »  Osx<2s 
p(x) = 
0 ; 2s<x<1 
(b) for 1/2<s<1, 
: 0 » Osxs(2s-1) , 
p(x) = 


C,/2(1-s) ; (2s-1)<x<1 ‘ 


These equations specify a series of nose-loaded or tail-loaded profiles, 
depending upon the value of s selected, as illustrated in the diagrams of 
Figure 4. When s=1/2, the profile is uniformly loaded over the entire 
chord. The contribution of the prescribed pressure distribution to the 
design lift coefficient is mC, This form of the prescribed pressure 
distribution also illustrates the general nature of the function p(x) 


which can always be written as 


p(x) = CG, fC) 
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Thus we see that the integrals have integrands which depend on the 


"normalized" pressure, f(x), which is defined to satisfy the condition 
f(x) dx = 1 


As a result, the step of factoring the quantity Cc from Equations (55) and 
(56), or more generally in any of the equations of the section on "The 
Sequence of Calculations," leads to no difficulties. 

When rectangular pressure distributions are used in Equations (55) 


and (56), they reduce to 


y,@) 


Cc 
= x, + (Vx -x,) eee = { (1-m) P(x) + mg (x,»8) 


ak. 1472s 
D(m,s) = (1-m) 272 = 2n(+v2) , m 1-/28 

ui T 
n ityv2s-1 


1-/2s-1 


a Q 
1+/2s-1 


F(x.) = ((42V7K)) Vote. - x [372 - an(14v2)) - tn(xt!4 4 Viws)) 
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Moreover for O<ss<1/2 we have 


2s-x 
g(x,.s) = ee Zi i < gp ees osm L gn (Ve) = ® anda + VIs/x,) 


=e 2s 

a £n alg > in ae 
° o 

and for 1/2<x<1 we have 


x 
“ fo} 1 1 
B(x,»8) Tighe ee eee 


2s-1 
- 1) + £n (14+V2s-1) + £n(2s-1) - 
Liv2s-1 Vx, aires 


1 1 1 2s-1 = 
age or eee * tx Ges) “™ x, : * wey tS Ay 
oO 


ea: See 
2(1-s) 


These closed-form expressions can be evaluated in accordance with the 

discussion following Equations (55) and (56). In order to determine a 

range of possible uw values for prescribed values of Cc. /T, s and Xo» it is 

sufficient for us to evaluate Equations (57) and (58) for m=0 and m=1. In 

the course of these evaluations, we must be certain that y, (1)20 is satisfied. 
In this context, we note that Fquation (57) shows for a fixed positive 

value of D, corresponding to prescribed values of s and m, that y,(1)/T 

decreases as Cc, /T increases. Thus the greatest allowable value of c,/T 

is limited by the condition y,(1)/T=0 in accordance with Equation (54). These 

limiting values of c,/T can be calculated for a range of m and s by finding 

the reciprocal of D(s,m). Results of such calculations are plotted in Figure 5. 


It may be recalled that when y,(1)=0 and m<1, the second design procedure 
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is being employed and that in this circumstance, occasions were found [4] for 
which no solution could be obtained. This situation was attributed to the 
fact that Cy and T could not be prescribed with complete independence. The 
curves of Figure 5 show more quantitatively the nature of the limiting relation- 
ship between Cc, and T which had only been described qualitatively in Reference 
4. Moreover, when m=1, the shockless entry condition prevails and a limiting 
form of the first design procedure applies because y, (1)=0 also. The condition 
y,(1)*0 shows that the first design procedure also has a limit for permissible 
values of the ratio Cc, /T. The range of parameters used in the first procedure 
did not exceed Cc, /T=.15, in Reference 3. Evidently, this range was not 
broad enough to encounter the limit illustrated in Figure 5. Further, Figure 5 
shows that the first design procedure permits a broader range of C,/T values 
than can be used for the second design procedure. We notice also that the 
influence of peak pressure location is greatest in the case when m=l. When 
m=0, the values of s have no effect on the permissible range of Cc, /T. However, 
aside from the case of m=0, we see that the designer will have more 
latitude with respect to Cc, /T values for tail-loaded profiles than for nose- 
loaded hydrofoils. In any event, the design values of C,/T must be selected 
to be less than those shown in Figure 5 for the various rectangular pressure 


distributions. In terms of the third design method, values of c,/T must 


be in the region of Cc, /T vs s of Figure 5 which is bounded by the curve 


m=] anc the C/T axis. For then, as we see from Equation (57), O<y (1) /Ts1/2 


and Equations (54) will certainly be satisfied. We note also that the m=1 


curve has an asymptote at x=.659 indicating for tail-loaded profiles that 
we can often choose Cc, /T to be as large as we please provided we are willing 


to move in the direction of shockless entry (m=1). Notice also that these 


results are completely independent of the values of x, or u. 


ee 


“~ 
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Having surveyed various limitations on the specification of the design 
point, we now turn to the influence of these factors on permissible values 
of uw. As noted above, for prescribed values of c./T, s and Xo? ic will 
suffice if we consider the two limiting cases m=0 and m=1 provided only 
that y,(1)20. First let us consider the case m=0 which corresponds to the 


flat plate profile with a point-drag singularity at an attack angle, 
a= y,@) + 2c, /1 . 


as can be seen from Equation (46). Clearly, the values of UU are independent 
of any prescribed pressure distribution and therefore of the parameter s. 
The results of calculations from Equations (57) and (58) for three values 
of x, are plotted in Figure 6. The locus of points for which y ,(1)=0 

is shown as the dashed vertical line in this figure. Only values of C/T 
to the left of this line are permissible. 

When m=1, Equations (57) and (58) can be evaluated again with the 
results plotted as shown in Figure 7. The three values of x selected for 
these calculations are the same as those used in Figure 6. For each of these 
Xo values, curves for nine values of s are shown with the corresponding boundary 
for y,(L)=0 shown as a dotted line intersecting the constant -s curves. For 


a prescribed value of s, only points to the left of the intersection give 


permissible yw and c,/T values. Figure 7 also shows that the lines from each 
value of s converge to a common point on the wU-axis. The location of this 
point on the » axis depends only on a This ordinate is determined from 


Equations (54) and (55) when C, /T=0. Thus, 


1 * (x, + vx )/2 . 
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These focal points are also plotted on the uw axis in Figure 7. 

The preceding trends, obtained from rectangular pressure distributions, 
have been considered chiefly because the analysis in these cases is greatly 
simplified and the exploration of basic theoretical trends and limitations 
is facilitated. On the other hand, comparison of these results at zero 
cavitation number with corresponding results for semi-elliptic pressure 
distributions [4] gives some indication of the effect of pressure-distribution 
shape on corresponding values of ut. This prescribed pressure distribution 
is illustrated for nose and tail loaded profiles in Figure 8 and its 


analytical description is 


(h/s)/x(2s-x) 9 h=2C, /ts for O<s<1/2 
p(x) = 


[h/(1-s)/(Q-x) (1-2stx) =, h=2C, /w(1-s) for 1/2<s<1 


Results obtained using this distribution are consistent with those of 
Reference 4. 

Comparison of various quantities in Figure 8 with those of Figure 
4 shows that the nomenclature for the rectangular and semielliptical pressure 
distributions has been made mutually consistent in order to facilitate 
comparisons of results from rectangular and semielliptical shapes. 
Figure 9 provides such a comparison at K=0 for values of the other parameters 


which are most sensitive to the details of p(x). These results were 


obtained numerically from Equations (55), (56) and (59) and from Equations 


(57) and (58). The curves labeled C, /T=2 and m=l1 show the variation to 
be expected in ) resulting from these two pressure distributions when 
x mel. For example, when s=.2 these curves show a difference in of about 


10%. For more rearward locations of the peak pressure, the difference in 
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the two values of pp is seen to decrease and to change sign for tail 


- 


loaded sections. The line for C,/T=0 shows U values which are independent 


of pressure distribution shape. Thus we expect variations in wy caused 


$e 


by shaping of the prescribed pressure to decrease with decreasing C/T 
and with decreasing m. 

A second comparison of the effect of pressure distribution shape 
is shown by Figure 10. In this figure, limiting values of Cc, /T are 
compared for semi-elliptic and rectangular pressure distributions. Again, 
only the value of m=1 is considered in order to emphasize differences 
which arise due to the shaping of otherwise similar pressure distributions. 
The curve labeled "rectangular" is simply the curve from Figure 5 for m=1l. 
The curve labeled "elliptic" was obtained numerically from Equations (47) 
and (59) with m=1 and y ,(1)=0. The asymptote for the rectangular distribution || 
is at s=.659. The asymptote for the semi-elliptic case was not determined. 
However, we believe that the asymptotic value of s lies slightly to the 
right of the value for the rectangular case because the: two curves cross 
at a high value of C/T which is off the graph. We also believe that the 
asymptotic value of s will be less than s=.7. The important point illustrated 
by Figure 10 is that the shape of p(x) affects the permissible Cc, /T range 
within which the constraints of Equations (54) can be satisfied. Taken 
together, Figures (9) and (10) suggest that there are two factors which one 
needs to consider when laying down profile design specifications. The more 
important of these is the chordwise center of pressure location for the 
prescribed distribution. We take this view because at m=1.0, shockless 
entry exists and the value of s is the center of pressure location for both 


rectangular and semi-elliptic shapes. The less important factor is that details 


of the pressure distribution shape will also affect the range of parameters 


pees 


ee 
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for which one can expect to obtain satisfactory profile designs. But 
this effect is seen to be considerably less pronounced than the primary 
influence of center of pressure location although it can lead to some 
uncertainty in the initial formulation of design specifications. 
Data for Semi-Elliptic Pressure Distributions 

The foregoing considerations show that the designer must accept certain 
constraints in the specification of the design parameters c./T, UH, Xo? K 
and p(x). Moreover, he can not be expected to arrive at a consistant set 
of specifications without some knowledge of the way these factors interact. 
In order to arm him with such knowledge, we present in this subsection a 
series of graphs which apply strictly to the semi-elliptic pressure distri- 
bution but which can also be used as a qualitative guide for other 
distributions. Unfortunately it will be expedient for us to restrict 
ourselves to the case of zero cavitation number. However, since foil 
design at K=0 is of considerable practical interest, these limited results 
are worth presenting. The reasons for restricting these data to zero 
cavitation number can be seen from the form of the general equations when 
the cavitation number is greater than zero. 

For example, we can start the general calculations for consistant 
parameters with the formula for the determination of the limiting values 
of Cc, /T. This result is obtained from the condition that y(1)=T as 
expressed by Equation (31) with the added condition that y,(1)=0. Equation 
(30), which corresponds to the condition n(1)=0, is used to eliminate a 
from Equation (31). In order to show how the term y.@) enters the 


resulting equation, we replace Equation (31), with Equation (30) used to 


eliminate a, and write 
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ativan 


-11 et+ad 
(1-m) ape F(a) + (a - 2 tan = 


22(6-ae) 


1 1 
Gr2z)tG) 4. 4 ben at tant 45 | f (x) dx 


ethoxy 2 5 VECRx) 


1 1 a 
“s | f (x) Fi ses dx - | £(x)&n sixty ix dx 
0 = 0 avx-V2-x 


Once the cavity length parameter, a, has been determined and y,Q) set 
equal to zero, the limiting value of Cc, /T can be determined from Equation 
(60) for values of m in the range, O<m<l. Note again that the pressure 


distribution p(x) can be written in the form 


p(x) = Ch f (x) 


where the function f is required to satisfy 


1 
| f£(x)dx = 1 
0 


The last integral is made equal to unity simply by scaling all ordinates 
of f uniformly if the integral should be found to differ from one. There- 
fore the three integrals in (60) depend only on f(x), the shape of p(x), 
and not explicitly upon the value of Cr because this factor can now be 
factored out of the right hand side of Equation (60) as shown. 

Before we can caluclate the limiting value of Cy» we must find the 
cavity length parameter a for prescribed values of m and K. This can 


be done if we put ak=y (1)=0 in Equation (40). For then we have 


A et I ET 
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1 
Cc 
1 L (2a* +3)e+ad , m Rx 
Bel = * oy 1 (1-) SS TH | oe f (x) dx A (40a) 
21 (1+K) 2(6-ae) R | x 
¥14+K 0 


This equation can be used to eliminate C, from Equation (60) after we put 


L 
y, (1)=0. Then we have 
i 
} 8 2a va “iva 1, _etad__ 7 (2-2x) £ (x) 
+ 
1 
; "* aha +2 tan > ay eee a 1] as tx & 
R 2 a | Gea) . va 
F 0 
t 1 etn tes 1 
F V2-x + avx . (2a 34.3) e406 4m{- (Eos 
- | f(x) &n Vix. ale dx| p=(1-n) ea a +7 | * £(x)dx . (61) 
0 0 


Since m, B(K) and T are specified, the cavity length parameter a is the 
, only unknown in Equation (61). Equation (61) can now be solved for the 

parameter a by iteration. Once values of a have been found for m=0 and 
m=1, say, Equation (60) can be used to determine limiting values of 
c,/T, corresponding to these two limiting values of m. 
Having found the limiting values of C/T for prescribed values of m 
and 6/T, one can now turn to the calculation of U for various values of 
m, 8/T and for values of C/T less than or equal to the limiting values. 
The first task requires that we determine the value of a when m, 8/T and 
C/T are specified and the normalized pressure distribution f(x) is given. 
This determination requires that we solve Equation (36) by trial and error. 
Once appropriate values of a have been found, Equation (32) can be used 
to solve directly for wu. 


By proceeding in the manner described above, one can calculate limiting 


values of Cy and consistant values of » for various values of C,, T, K, m 


gy 


Far 


— i theeeene! 
. 
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and Xo" Having the limiting values of Cc. enables one to be certain that 
y, (1)20. However it is apparent from Equation (36) that although the 
quantities c,/T, B/T and m are involved, the quantity reer ce ii also 
appears. Evidently it is not possible to reduce the problem to one 
involving Cm 8/T and wu as primary variables except for the special 
case of K=8=0. In this special case we have already seen that C/T and 
yu serve as useful parameters for the correlation of consistant values of 
m, X, and characteristics of the pressure distribution shape. Therefore 
we will give further examples of permissible ranges of yp for K=0 because 
this is a practically important case and because the trends displayed 
will have qualitative value for the interpretation of results at other 
values of K. 

For the case of K=0, we can use the special formulae of Equations (55) 
and (56), with Equation (59) for p(x), to calculate the set of consistant 
values which satisfy Equations (54). The first set of results is shown 
in Figure 11 which gives the limiting values of C/T for m=1 and m=0 for 
values of s between .1 and .5. These values of s are more than sufficient 
to cover the practically useful range of c,/T. Figure 12 shows values of 
u for y,(1)20 for three values of Xe when m=0 and K=0. Figures 13, 14 and 
15 show consistant values of ww when K=0 and m=1 for the three values of Xo 
given in Figure 12. In these three figures, it is also necessary to show 
the effect of the peak pressure location s. As we have seen already, there 
is no dependence on s when m=0. In all of these figures the values of 
y,(1)=0 for each Cc, /T correspond to those shown in Figure 11. 


It is found by comparing Figure 12 data with those from Figures 13, 


14 and 15 that there is only one permissible value of w for all values of 
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m in the range O<msl] when s is tees to .2. If s is less than .2 the value 
of yt at m=0 is greater than that at m=]. On the other hand, if s exceeds 
-2, the value of » at m=0 will be less than that for m=l. For any 
particular design-value of C/T the permissible range of wu will lie between 
the value for m=0 and the value determined by m=1 and the design value of 
s. We also note that the cut-off values of Cc, /T for y,(1)=0 differ when 
m=0 or m=1. This is illustrated in Figure 11 for the semi-celliptic pressure 
distribution and in Figure 5 for rectangular distributions. This fact means 
that these preliminary data for semi-elliptic distributions are not entirely 
complete even for the zero-cavitation number case. There will be cases 
in which the design value of Cc, /T will require that the lowest permissible 
value for m (=m', say) must be greater than zero. This simply means that 
the permissible values of m will now lie in the inteval 0O<m'<m<l if y,(1)20. 
In spite of the uncertainty remaining in Figures 13, 14 and 15 regarding 
the value of corresponding to the lowest m-value for the design value 
of s, we can still start the calculation by the third design method using 
the data from Figures 12, 13, 14 and 15. If one finds that y,Q1)<0 for 
some value of wu which was thought to be permissible this case can be 
discarded in favor of a value of Lt nearer to that for m=1. 


Alternatively, one can assume that the design value of C, is the 


L 
limiting value. Then, when y,(1)=0, Equation (55) shows that the 


necessary value of m can be determined by 


372 - &n(1+/2) - nT/C, 


3V2 - &n(1+v2) + =. 5 | E(x) on 2A ax 


9 1-/7x 
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We can then put m=m' in Equation (56), along with prescribed values of 


C,, T and Xo in order to find the corresponding limiting value of y. This 


1 Ph 
value of wu and the value obtained for m=1 determine the admissible range 


of uw for the case when the design value of Cr lies between the two 


limiting values of Cy for m>m'>0 and m=1. 


The Design Procedure for Arbitrary Cavitaton Number 


For nonzero cavitation numbers, it is not possible to present the kind 
of preliminary data as compactly as we did for K=0. The primary reason 
is that now both B(K) and the quantity (1+K) appear in the equations. 


Therefore these input parameters together with other prescribed quantities 


which are C_, x_ and T, in addition to the pressure function p(x), seem 
L ° 


to make the task of reducing the number of governing parameters not too 
useful. Moreover, the selection of values for the several input quantities 


must always be such that the constraints, 


y.@) > 0 


are satisfied. For otherwise, our calculations could represent a physically 
unrealistic solution, as we have noted previously. 


We have already seen in the case of K=0 when C_ and T are prescribed 


L 
that the constraints of Equation (54) can lead to certain "cut-off" values 
for the parameter yw which correspond to m=0 and y ()=0. In these examples 
the shape of the prescribed pressure distribution is represented by the 


single parameter s. For nonzero values of K, we will continue to use this 


parameter for purposes of preliminary analysis. However, we will now 
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abandon the use of the parameter C/T and simply treat yp, C_, T, x and 


L? 
K as input quantities along with the value of s. Now we will envisage 
fa the constraints defined by Equations (54) as providing an admissible 
region ina Cy-u plane. The extent of this region will depend on the 
values assigned to K, T, X and arbitrarily, it will encompass values of 
s between .1 and .9. The admissible region of the Cu plane is 


determined by Equations (30), (31) and (32) which follow from 


n(1) 
y (1) 


y(x,) 


respectively. The boundaries of this region are found from the constraints 
of Equation (54). 


If one prescribes a design value of Cy in addition to values of K, 


Es Xo and s three different cases can be represented in the C,-y plane. 


L 


These cases, which are illustrated schematically in Figure 16, arise from 


the fact that the limiting values of C for m=0 and m=] will be different. 


Such differences for the special case of K=0 have been shown already in 


Figures 5 and 7. Similar behavior will also occur when K>0. As we 


have seen, the limiting value of Cy for m=0 does not depend on any 
prescribed pressure distribution. This is illustrated in Figure 12 
for K=O. Similar trends for Ut versus Cy occur when m=0 and K>OQ. For 
any cavitation number one can calculate the limiting value of Che 


regardless of the form chosen for p(x), once T, X, and K have been 


prescribed. This limiting value of Crs which corresponds to m=0 and 


y,,(1)=0, is denoted by i in Figure 16. In terms of the value of s 
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which characterizes the semi-elliptic pressure distribution, it is found 
in the C,-y plane that the slope of the ray for m=0 from the focal point 


L 


on the Cy axis is very close to the ray for s=.2 when m=1. Since s 


defines the abcissa of the centroid of the semi-elliptic pressure distri- 


bution one may suspect that other nose-loaded pressure distributions, 
having their centroids located at about 20% of the chord, will be similarly 


related with respect to the slope of the ray for m=0 in the C,-py plane. 


L 
When the ray for m=0 and the ray for m=1 coincide for a centroidal location 
s, a single value of yw is defined as long as the design Cy is less than 
or equal to CLo- This value of uw for such a value of Cy is denoted by 
Uo in Figure l6a. Physically this means that it is possible to prescribe 


pressure distributions p(x) and values of C,, K, T and X for which only 


L? 
one value of up is permissible regardless of the values of m and y,@) which 
can then be paired to produce a series of profiles of varying performance 
and geometry. This is a rather special circumstance which will not often 
be encountered. It is far more likely that the m=1 ray for the prescribed 
pressure distribution will lie above the m=0 ray or below it as illustrated 
in Figure 16a. For the permissible range of wu which can be prescribed 
when the design Cc lies well to the left of CL will then lie between Wy 
and Yo or Ug and Uo> depending on whether the prescribed value of s lies 
below the m=0 ray as is the case for S$,» or above the ray for m=0, as is 
the case for 8, in Figure l6a. 
With the foregoing qualitative description of the circumstances 

surrounding the case illustrated by Figure 16a in mind we now turn to the 
calculation required to implement the determination of permissible values 


of the parameter |i. The first task is to calculate the limiting value, 


Cho: This calculation can be carried out by means of Equations (60) and 
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(61) with m=0 and y,(1)=0. First, Equation (61) can be solved by iteration 
to determine the value of the cavity-length parameter, a, for the prescribed 
values of T and K. Then CLo can be found from Equation (60). Next when 
3 m=1, we can put y,(1)=0 and repeat the iteration of Equation (61) in 
order to find the value of a and then the corresponding value of Cc 


can be found from Equation (60). The value of Cy determined by this second 


A 
f sequence of calculations corresponds to the intersection of the ray Sy 
5 or Sy: with the solid curve labeled y,(1)=0 in Figure l6a. These limiting 
: C values are denoted by CLy and Cho in Figure l6a. 
| Now, in order to consider the case in which the design Cy is less than 
* the values CLo and Chy which have just been established, we note from 


Equations (34) and (34a) that the relationship immediately following Eq. (34a), 


m = 
f(a) = Fy + 714K) H, = 0 ; (62) 
does not depend upon x, OF Ue Therefore if design values of Cis K, T and 
p(x) are used in (62) along with m=0 and m=1 the value of the cavity length 
parameter can be found iteratively. Then the two pairs of quantities, 


(m,a), and the design quantities C_, K, T, p(x) and the remaining parameter 


ee 


xX,» can be used in the equation 


m 
T  ~ W(+K) Hy ae? 


in order to determine , and Wy or Mos The functions Ho and Fp are defined 


by Equations (35) and (35a) respectively. The permissible range of yu 


which one can specify for subsequent design calculations will then be in 


the interval (uy> Uo) or (gs Uy) as noted previously. 
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The foregoing discussion, while helpful for describing the basic 
nature of the calculation process for the estimating mode, indicates a 
rather cumbersome numerical procedure. It has been found more effective 


to program the estimating-mode calculations for Cho and CLy or CL to 


¥ follow a different path from that discussed above. A more efficient 

4 determination of these limiting values of C. can be based upon 

4 

ta Equations (62) and (63) because they are also needed for other parts of 


the design process. In particular, since Equation (62) does not 
depend upon X or wt, if m is chosen in advance this equation contains 
the unknowns C, and a. Then if we put E=y ,(1)=0 in Equation (40), 
thereby obtaining Equation (40a), we can solve it for Cc in terms of 
a, m, K and p(x). The iteration of a is then conducted by calculating 
an estimate for Cc for each estimate of a from Equation (62). This 
4 estimate for Cc along with the prescribed values of K, T, m and p(x) 
can be used in a regula-falsi ‘ routine to find that value of a which 
satisfies Equation (62). Clearly, this process will produce the limiting 
values Cho and CLy or CLo illustrated in Figure l6a. Corresponding values 
of Uo and Wy or Uy then follow from Equation (63). 

Of course, if the design pressure distribution is such that the ray 


for m=1 lies above the ray for m=0 as illustrated for s, in Figure l6a 


2 
the procedure just outlined for the estimating mode holds without 
modification for values of the design Cc less than or equal to Cho: However, 


suppose that the design value of C lies between Ch and CLo and that 


the value of s is such that the ray for m=1 is below that for m=0, as 


shown by the ray for s, in Figure 16a. Then the calculations for the 


1 


estimating mode must be modified in accordance with the illustration of 


Figure 16b. In this case the value of Uo can be found from the procedure 
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described above when m=0 and the design values of Ce K, T and x, are 
prescribed as usual. However when m=1 the design value of C exceeds 
CL so that this value of m is not available for profile design. Instead, 
we must find a value of m which is less than unity but which is just large 
enough to permit the design C to be the limiting value of C. for y,(1)=0. 
This value of m=m' will lie on a ray between m=m‘' and m=0 as shown in 


Figure 16b. The value of m' can be found from Equation (40a) by solving 


it for m. Thus we find 


2mB (14K) /C, - [ (2a7+3) e+a6]/2(S-ae) 


1 
+| EX ¢(x)ax - [ (2a7+3) e+ad] /2(8-ae) 
0 


Then since Cy in this equation is the design Ch. Equations (40b) and (62) can 
be used in a regula falsi routine to determine a and m'. Then m', a and 
the other design parameters can be used in Equation (63) to calculate 
u'. Permissible ) values will lie between y' and Uo: 
The third case for the estimation mode calculations arises when the 
design Cc lies between Cho and CLo and the pressure distribution leads to a ray, 
such as So in Figure l6c, which lies above the m0 ray. In this case, the 
value of Hy for m=1 can be found by iterating for a in Equation (62) and 
then solving for Uy from Equation (63). However, the value of m=0 is 
not accessible for profile design in this case. Instead we must find 
an intermediate value of m between 0 and 1. The lowest such value, m', 
which corresponds to the design value of Cc will be found when y,(1)=0. 
Thus we can use Equations (40b) and (62) in a regula falsi routine to 


determine limiting values of a and m', as noted above. Then the 
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corresponding parameter ]}' follows from Equation (63) and the permissible 
range of 1} lies between i’ and y+ 
It can be seen that, as the design Cc, increases and the other para- 
meters are held fixed, the value of s (or the centroid of the prescribed 
pressure distribution) must move nearer to the profile trailing edge. 


Indeed, if the design C_ is greater than CLo in Figures i6a and b, no 


L 
solution can be found. Of course, it is evident from calculations of 

the sort which lead to Figure 5 that this rearward motion of the centroid 
of the prescribed pressure does not continue indefinitely because it is 
found that a centroidal distance between .6 and .7 will lead to infinite 
values for the limiting Cy when m=1. For values of m less than unity, 
this asymptotic value moves farther to the rear of the profile. However, 


for practically interesting values of C, and the other design parameters 


L 
this restriction is not too limiting. 

In the design mode the calcula) ions can go forward using the permissible 
u values from the estimating mode. These calculations are carried out 


as described in the section on The Sequence of Calculations. In this 


case the iteration for a is executed using Equation (33a). 


THE THIRD DESIGN PROCEDURE -- SUMMARY OF RFSULTS 
Preliminary Remarks 
Nearly all of the calculated results to be discussed in this section 
are based on the semi-elliptical pressure distribution of Figure 8. Since 
the third design procedure contains the First and Second procedures as 
special cases, the numerical results presented here for the case of 


m=1, corresponding to shockless entry, will permit some regions on the 


wetted surface of the foil to be at the cavity pressure. This situation 


-58- 5 July 1977 
BRP: JF: jep 
violates the basic restriction that the pressure on the profile should 
be greater than the cavity pressure except at the trailing edge (and for 
shockless entry at the leading edge too.) Nonetheless, the basic boundary 
condition (vi) is satisfied. We could have avoided this contradiction if 
we had used the double~ellipse pressure distribution which was used for 
the first procedure in Reference 3. It was shown in that study that the 
hydrodynamic performance at ee oe entry resulting from the semi- i 
elliptical pressure distribution or the double-ellipse distribution 
exhibited only small differences. Therefore, we decided to overlook this 
contradiction for this one limiting case in the interest of simplicity 
and with the knowledge that valid design trends would still be obtained. 
All other cases, corresponding to m<l, involve the angle of attack or 
nose singularity so that the pressure on the profile will always exceed 
the cavity pressure except at the trailing edge. 
For the sake of completeness, all hydrodynamic performance trends 
calculated in the course of this study are tabulated in the Appendix. 
The numerical methods used and the computer program developed for the 
calculations are described by Mr. Fernandez in a separate report [1]. 
Hydrodynamic Performance Trends 
Perhaps the single most critical measure of hydrodynamic performance 
for fully cavitating profiles is the lift-to-drag ratio. Accordingly, 
our investigations have concentrated on the relationship between this (7 
factor and other design parameters. For example, Figure 17 shows the 
influence of peak pressure location on lift-to-drag ratio at zero cavitation 
number for a range of lift coefficients. For each lift coefficient, two 


curves are plotted corresponding to the largest and the smallest permissible 


value of j) when xml and T=.1. This illustration shows that the selection 
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of can have an important effect on hydrofoil performance depending 
upon the location of peak pressure. It appears that for s-values in 
the neighborhood of .3 the effect of uw on L/D is minimal and that for 
tail-loaded profiles , it is best to use the smallest permissible 


U-value, at least when K=0. For nose-loaded profiles the trend is 


reversed. However, it was found that maximum values of Cc, /T as 


illustrated in Figure 10 were exceeded at all Cy values except C,=.08. 
Thus no data are available at values of s equal to .2 and .1 for the 
higher lift coefficients. 

Figure 18 provides the same information as Figure 17 except that now 
the cavitation number is .2. The most important factor revealed by this 
illustration is the significant decrease in L/D-values when K is increased 
from 0 to .2. We also see that at lower values of Cc the effect of u 
on L/D increases for tail-loaded profiles as Cy is increased. For nose- 
loaded hydrofoil sections the effect of Cr changes is not large. 

Both Figures 17 and 18 show the consistant trend of increasing 
L/D with increasing Cc. This same trend was found for the first and 
the second design procedures [4]. 

Figure 19 shows the effect of peak pressure location on L/D when 
the cavitation number varies and Cy is fixed at .16. Again, T=.1 and 
xml This figure also shows that the effect of permissible ,) variations 
is not too great for nose-loaded profiles, but for tail-loaded sections 
it can be significant, especially at very low K values. 

Figure 20 shows the effect of increasing cavitation number on L/D. 
The value of s=.3 has been selected as a fixed parameter so that the 
variations due to the range of permissible !) values is small. From our 


preceding discussion we recognize that had we selected a tail-loaded 


a 
Me 
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profile instead of the particular loading given by s=.3, the effect of wu 
variations would have been much more pronounced. However, the chief 

trend illustrated is the decrease of L/D with increasing K. The present 
trend agrees with that found for the first and second design procedures [4]. 

The influence on L/D of cavity thickness, T, for the range of design 
C, values is shown in Figure 21. These curves are drawn for the fixed 
values of K=0, s=.3 and Xue As noted already the effect of variations 
in 4 is not great because of the s-value selected. The general trend 
shown in this graph is the dramatic decrease in L/D with increasing cavity 
thickness. This trend has also been observed for the first and second 
design procedures [4]. 

The chief new aspect of the present investigation is the introduction 
of the two-point specification of cavity thickness. For a prescribed 
value of T, the parameter 1) controls the specification of the cavity 
thickness near the nose. Further information on the effect of UU variations 
upon the values of L/D is given in Figure 22. In this figure the values 
of K=0, T=.1 and xed are held fixed. Corresponding values of L/D 


and yt are plotted at the two design C,'s of .08 and .16 for the range of s 


L 
values. Among other things, this illustration shows how L/D is 
independent of when s is somewhat greater than .3 but less than .4. 

We also see that in the neighborhood of s=.2 there will be only one 
permissible value of y) for the whole range of L/D within which m varies 
between 0 and 1. This figure also shows that in most cases L/D tends 

to decrease as ) increases. The exceptions to this situation occur when 
s lies between .2 and .3, roughly. We also see that at the higher Cc, 


value of .16 most of this region is denied us by the cut-off corresponding 


to y,(1)=0 and the permissible value of CL ax’ 
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Figure 23 is a cross plot of data from Figure 22 for constant values 
of wt. The purpose of this plot is to examine in some detail the L/D trends 
for various peak pressure locations for the two values of design Cy and 
for selected values of uu. The point at issue is whether or not one 
preserves the trends of increasing L/D as he prescribes peak pressure 
locations which move towards the nose when H is held fixed. As Figure 23 
shows, the answer to this question is yes in most cases, at least when 
K=0. On the other hand, for some C values this trend is reversed. This 
reversal occurs when the Ut value selected lies to the left of the 


vertex at m=0, as shown for C,=.08 in Figure 22. Indeed, as already 


L 
noted, if »p is selected so as to coincide with the vertex, one finds 

only one permissible value of s for the entire range of L/D corresponding 
to O<m<1. In this case the trend line in Figure 23 would be vertical. 
Thus, the idea that nose loading is best from a purely hydrodynamic view- 
point is not always valid as already noted in Reference [4]. However, 

it appears to be true in most cases so that we might regard it as a 

useful rule of thumb. Another point illustrated again by Figure 23 is 

the fact that for constant UL values, the permissible range of s is limited 
by the constraints on m and y,@). 

Figure 24 is similar to Figure 22 except that now K=.2. The effect 
of this increase in cavitation number on the trends described above is 
remarkable. In many ways the situation is reversed. This is definitely 
true with respect to the effect m on L/D for nose-loaded profiles. For 
example, when s=.1, L/D increases as | increases and at the same time 
m decreases from 1 to zero. We do see that near s=.2 there is only one 
value of w for all values of L/D corresponding to m values between 0 


and 1. In both plots of Figure 24 m is zero at the vertcx where all 
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rays intersect. Then as a point moves from the vertex to the end of a 

4 ray, the value of m increases from 0 to 1. The plot for C,=-08 and K=.2 
fi is shown to a larger scale in Figure 25 so that the effects of y) upon L/D 
can be seen more clearly. Evidently, at this value of cavitation number . 
Fe and at this low value of Che better L/D's can be achieved if s is between 
-6 and .7 and m=1. On the other hand if c,=.16, the best value 


+e of s is near .4 with m=l. However, even then, the L/D values obtained 


are low. 
In the course of our numerical studies we carried. out some computations 
ae 
i to test the sensitivity of L/D to the accuracy of the lowest and highest 


permissible values of wu. It was found for x et that if L/D and m were to 
be determined to four significant figures, these limiting values of uw needed 
to be known to five significant figures. 

The preceding data have been presented in order to highlight some of 
the more important performance trends exhibited by third-design-procedure 
data. In order to present a somewhat more systematic overview of hydro- 
dynamic performance, we present in Figures 26 through 35 plots of L/D 
versus Cy: These performance maps show contours of design Cy and s for 
cavity thickness values of T=.10, .15 and .20. Separate maps have been 
prepared for five values of the cavitation number, K=.00, .05, .10, .15, 
-20. For each value of K individual maps are given for the greatest and 
the least values of , consistant with the constraints O<ms<l and y,(1)20. 
Since the value of the greatest yp and the least » will vary widely from 
point to point on each map, the two maps for each K-value are simply 


designated by "MU MAX" or "MU MIN" on the plots. All maps have been 


contructed for Xcel. Complete tabulations of hydrodynamic data for 


these ten maps are given in the Appendix. Therefore one can look up other 
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pertinent data not presented on these plots should he need more detailed 


information. These maps seem to be rather self explanatory. Therefore 
we will only comment that they show again that nose-loaded profiles most 
often seem to have larger L/D values than tail-loaded sections. They 
also show ranges of design parameters for which the third design procedure 
can give a solution. Of course, not all cases investigated are presented 
in these performance maps. However, in order to present the hydrodynamic 
data behind all plots presented in this study these additional results 

are also given in supplementary tables in the Appendix. 

In Reference 4, it was found useful to introduce the concept of 
“hydrodynamic equivalence". In this concept the design point (K, Ch. T) 
and the center of pressure location, x, are the primary factors which 
govern hydrodynamic performance. It was found that if these parameters 
were the same for two profiles, their performance would be closely the 
same even though the shape of their pressure distributions is different. 
This concept was found to be useful even for comparing the performance 
of profiles designed by the first and second design procedures. Therefore 
we would expect this concept to be of value for comparisons of profiles 
designed in accordance with the third design procedure. In fact the 
center of pressure location is given in all tabulations in the Appendix 
so that one can test the concept of hydrodynamic equivalence again. We 
have not made this test in the present study. Instead we have preferred 
to explore the consequences of the third design procedure on the profile 
geometry. As we found in Reference 3, profile geometry is sensitive to 
the shape of the pressure distribution while hydrodynamic performance is 


less so. 
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Profile Geometric Trends 
Typical profiles resulting from the third design procedure are 
illustrated starting with Figure 36. In all of these illustrations the 
upper graph shows to scale the profile wetted surface and the upper cavity 
contour in accordance with the nomenclature defined in Figure 1. In 
these figures the chord lines of all profiles are on the x-axis from 


x=0 to x=1.0. The y-axis shows contour ordinates to the same scale as 


the abscissa of each point. Note that only that portion of the cavity 


upper surface which lies directly above the wetted surface of the profile 
is plotted in these graphs. The trailing edge of the wetted surface is 
located at the point x=1.0 on all plots, and the lower surface of the 
cavity, which would extend beyond the trailing edge is not shown at all. 
The cavity thickness, T, for each design can be measured directly from 
the graphs as the vertical distance from the trailing edge of the profile 
(1.0, 0) to the upper surface of the cavity (1.0, T). The space between 
these two contours from the profile nose at (0,0) and on into the cavity 
beyond the trailing edge of the wetted surface if need be, is available 
for hydrofoil structure. We have not drawn in any foil structure in 
these plots, preferring to restrict our attention to the hydrodynamics 

of these flows. Of course we would assume that there would be some 
clearance between the upper surface of the cavity and the uppermost 

parts of the profile structure at the design condition. The lower plot 
in each figure shows the pressure distribution along the chord of the 
wetted surface with positive pressures plotted downward. These plots 
show the net perturbation pressure due to the prescribed pressure distri- 


bution plus the "flat plate" distribution due to the angle of attack above 
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the ideal value for shockless entry, as discussed in the text following 
Equation (43). 

The examples selected for study have not been chosen because the 
profiles exhibited have good or poor performance or show favorable or 
wibayoeatla geometry. Instead our purpose has been to illustrate 
typical geometric trends which follow from changes in key hydrodynamic 
parameters. In this regard, we note that many of these general trends 
have already been discussed in Reference 4, where profiles designed by 
the first and second procedures are illustrated. Therefore the present 
discussion centers on the chief new ingredient contained in the third 
design procedure, namely; the influence of u, the cavity thickness para- 
meter at the nose, upon geometry and net pressure distribution. 

For example, Figures 36 and 37 show two designs at zero cavitation 
number having the peak pressure of the prescribed semi-clliptic pressure 
distribution at s=.2. Figure 36 corresponds to the limiting case of m=l 
so that it is a "first-design-procedure" hydrofoil section. We have 
noted previously that the value of s=.2 corresponds closely to a case as 
illustrated in Figures 22 and 23 for which » is near the vertex corresponding 
to m=0 in a plot of L/D vs u. Accordingly the value of wu changes very 
little over the permissible range of m. In this case, uw values shown 
in Figures 36 and 37 differ only in the fourth decimal place. The value 
of y,Q) for the profile of Figure 36 is .0081 so that the nose in this 
design is almost sharp. The value of y, (1) for the foil of Figure 37 
is zero and the nose is sharp. The value of m is .6382 for this case so 
that the pressure distribution shows significant flat plate contribution. 


Because of this fact, the reduced amplitude of the prescribed semi- 


elliptical pressure distribution leads to less camber of the wetted 


| 
i 
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surface when compared to the case for m=] shown in lee 36. The 
corresponding difference in lift-to-drag ratio for the two cases is 
seen to be fairly small. Of course, neither of these designs is 
suitable for practical application. There is little or no clearance 
between the wetted surface and the upper-cavity contour near the profile 
noses. In addition, the design of Figure 36 shows an extensive region 
of zero pressure on its wetted surface. This situation is easily 
remedied by use of an additional quarter-elliptical distribution of low 
intensity distributed over the interval which now has zero pressure, as 
was done for the first procedure foils of Reference 1. The reason that 
this measure was not adopted in the present study is because our interest 
centers upon designs in which the flat plate contribution will cause the 
eke surface pressure to exceed zero everywhere between the leading and 
trailing edges. The limiting case for m=] is already treated adequately 
in Reference 1. The resulting comparisons of performance and geometric 
trends are influenced very little by this minor violation of one of our 
basic constraints. Of course one can not recommend such a prescribed 
pressure distribution for practical applications. 

The next five illustrations, presented in Figures 38 through 42, 
illustrate the case in which s=.3 and K=0 and one expects to see a 
greater permissible range of u values. We expect this because as can 
be seen from Figure 22, when s=.3, a considerable range of ) can take 
place with relatively small changes in L/D resulting. In fact if one 
can select an s-value which is on a horizontal line through the vertex 
at m=0, then L/D wll not vary over the permissible range of yp. In 


Figures 38 through 42, we have 19.906<L/D<20.135 so that the condition 


of L/D invariance with respect to y is nearly satisfied for the 
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permissible range, .1464sys.1577. The profile and pressure distribution 
shapes show descernible variations in this range. These variations can 
be understood quite readily if one examines the values of m and y,Q) 
for each of these profiles as given in Table I. For then, the general 


discussions 


Table I 


Fig. 38 39 40 41 42 
m - 9998 - 8537 «1075 - 5609 -4148 
y,@) - 0200 -0150 - 0100 -0050 - 0000 


which relate to the importance of these factors on profile and cavity 
shapes as noted for hydrofoil sections designed at s=.2 are applicable 
to this case also. 

Continuing our study of geometric trends at K=0 and T=.10 in Figures 
43, 44 and 45, we show profiles designed for s=.7. For such tail-loaded 
sections the values of y,@) needed to provide the required cavity thickness 
are higher than is the case for nose-loaded designs. Therefore, these 
profiles show higher drag and more nose rounding. As expected, the region 
of under camber moves to the rear and the moment coefficients of these 
profiles are greater in magnitude than is the case for nose-loaded 
designs. The values of m and y.@) for these profiles are given in 


Table II below. 


Table IL 


43 44 45 
- 9998 «6079 -2158 
- 0516 +0258 -0001 
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In Figures 46, 47, and 48 we return to nose-loading with s=.3. 
The values of wt} selected for these graphs have been chosen to 
produce values of m as close as possible to .5 in order to retain 
pressure distributions which are the same for the three cases. The 
values of m were .4995, .4994 and .4994. Therefore the changes of 
profile shape are due to input parameters other than the pressure 
distribution. That is, we have made certain in these comparisons 
that x, the center of pressure location, does not change. For example, 
Figure 46 shows the effect of changing cavity thicness for K=0. The 
case shown has an m-value between those of Figures 41 and 42 as indicated 
in Table I. An interpolation based upon the values in Table I indicates 
that y,(1)=.0029 for m=.4995. These values pertain to T=.1. In Figure 
46, T=.2 and y,(1)=.0529, which is about 18 times larger than the value 


of y,@) when T=.1. Thus the change of cavity thickness from .1 to .2 


ean be ascribed entirely to the increase in the value of y, (1). As a 


result of this increase in y, (1) there is rounding of the profile nose. 
In Figures 47 and 48, we return the cavity thickness, T, to its 
previous value of .1. However, the value of K is increased to .1 and 
then to .2. Again, because m is held constant, the cavity thickness 
requirement of T=.1 for these K values results in increased values of 


y,@) as shown in Table III. 


Table III 


47 
0 el 
- 0029 -0190 
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Again the increased rounding of the profile nose which accompanies this 


increase in y, 1) is observed. 


As our final study of profile geometry to be expected from designs 
developed from the third design procedure, we consider some sections 
based upon a modified "three-term" pressure distribution [4,6]. This 
distribution is shown in Figure 49 for a design lift coefficient of .12. 
As illustrated in the figure, two uses are made of this distribution. 
First we reverse the distribution so as to produce a nose-loaded profile. 
Then we use the distribution in the customary way to provide a tail- 
loaded section. 

Four profiles, resulting from the third design procedure when 
T=.1 and K=0, are shown in Figures 50, 51, 52 and 53 for reversed-three 
term and three-term pressure distributions. We believe that the study of 
hydrodynamic equivalence reported in Reference 4 supports the view that 
the value of peak pressure location, s, recorded on these plots is no 
longer a significant parametric quantity for comparison of performance 
with results from semi-elliptic pressure distributions, for example. 
Instead the important parameter as far as hydrodynamic performance is 
concerned is the center of pressure location, x. These values are 
tabulated in the Appendix. On the other hand, since profile geometry 
is sensitive to details of the prescribed pressure distribution and the 
point-drag amplitude these distributions tend to be in a class by 
themselves with respect to the geometry resulting from the application 
of the third design procedure. We might consider these profiles as 
being more representative of the kind of pressure distribution one 


might wish to specify for actual use. 
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Figures 50 and 51 show the profiles derived from the third design 
procedure for the reversed three-term prescribed pressure distribution. 
In Figure 50, m=1.0 and y, (1)=.0228. Figure 51 shows the profile and 
net pressure distribution obtained when m=.3840 and y,(1)=0. Therefore 
the value of uw shown on Figures 50 and 51 are —_— and a respectively. 


Note that the L/D value for the profile of Figure 51 exceeds slightly 


& that shown in Figure 50. This increase is caused by the fact that the 

5 location of x moves forward because of the flat plate solution, shown 

Ps to be part of the pressure distribution of Figure 5l. 

* Figures 52 and 53 show the geometric trends obtained for the three- 
term prescribed pressure, corresponding to tail-loading. Again, Figure 

n 52 corresponds to the limiting case for m=l and Figure 53 to the limit 

| for y,,(1)=0. In this latter case, the value of mis .2274 so that 
q x is moved towards the nose compared to the profile of Figure 52. 


Accordingly, the value of L/D is increased for this case as compared 

to that shown in Figure 52. Again the values of — shown for these two 
profiles are limiting permissible values, Ul ax’ 2085 and ul in 2543. 

It is worth noting that the value of L/D for the profile of Figure 53 

is nearly as good as that for the shockless entry case of Figure 50. These 
two cases show again the important influence on performance of the center 
of pressure location as distinct from peak pressure location. On the 


other hand, the effect of distributing the prescribed pressure 


gradually over the chord instead prescribing the more concentrated peaks 
as obtained in many of the semi-elliptical examples given previously, is 
shown by the more gradual camber of the wetted surfaces in these last 
four examples. Perhaps, one conclusion which can be drawn from these 


comparisons is that the designer can, by a judicious selection of 


+e 


= Tile 5 July 1977 
BRP: JF:jep 


prescribed pressure distribution and the value of yt, obtain profiles 


of better than average performance and still provide sufficient clearance 
between the upper surface of the cavity and the wetted surface to allow 
realistic structural criteria to be satisfied. The perusal of this 
objective is beyond the scope of the present study. 
Off-Design Performance | 
Once a profile has been designed, it is important to examine its 
performance at possible operating conditions which differ from those 
considered in its design. Considerations of this sort depend upon the 
solution of the direct problem and the literature is replete with various 
approaches to its solution. One approach which appears to work well with 
the present inverse technique is that described in Reference 4. This 
particular linearized approach has been designed to permit one to 
determine the operating range of a particular design in terms of cavitation 
number and minimum attack angle for cavity profile interference. 
In order to illustrate the off-design performance, one might expect 


to achieve with the third design procedure we have selected one profile 


from those presented in the preceding section. The profile of Figure 40 


has been analyzed for off-design performance and the results are 


presented in Figure 54. In this figure, the lift-to-drag ratio is plotted 
against the cavitation number for selected values of attack angle, a, 
and cavity length, &. Regions of cavity-foil interference have also 
been delineated on this figure. In arriving at the boundary between 


interference and noninterference, we used a rather simple shape for 


the upper surface of the hydrofoil section. It was assumed that, t, 
the ordinates of the upper surface of the hydrofoil is 6 times the 


ordinates of the upper surface of the cavity at the design point. ‘hat 


ay i et 
t 
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is, t=O0y(x). In the example of Figure 54, 6 is 80% and the hydrofoil ) 
has a kind of wedge shape. The boundary marked in Figure 54 defines 
the locus of points at which the upper surface of the cavity at off- i 
design is just tangent to the upper surface of the foil at some 


chordwise point. 


CONCLUSIONS 


The work completed to date has provided a rather flexible 
preliminary design tool which one can use to explore the important 
hydrodynamic parameters surrounding the design of a supercavitating 
hydrofoil section. The work has been carried out with the needs of 
the designer who wants to explore mixed-foil concepts in mind. Naturally, 
we have limited our considerations to the hydrodynamic aspects of the 
design process. We hope that the methods proposed will enable the 
important compromises between hydrodynamic and structural criteria to 
be addressed more efficiently than might otherwise be possible. 

It has been found that the results of the present study reinforce 
conclusions of our work on the first and second design procedures with 
respect to the effect of center of pressure location and other parameters 
on hydrodynamic performance and the influence of the prescribed pressure 
distribution shape on wetted surface geometry. 

The chief new result from the third method of profile design results 
from the two-point cavity thickness control which is possible with this 
more general theory. In calculations based upon this theory, one 
specifies the cavity thickness at a point Xo Near the nose to be uT, 


where <1. Because of the fact that the pressure on the wetted surface 


must exceed the cavity pressure, it is found that » can not be prescribed 


om “e 
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with complete arbitrariness. The permissible range of yw is found to 


depend on T, x? C_, K and, to some extent, upon the shape of the 


L? 
prescribed pressure distribution. In spite of these restrictions on the 
range of wu, the addition of this parameter to those already available 
from the first and second design procedures significantly enhances the 
flexibility of the third design procedure. 

Of course this procedure, being based upon linearized theory, like 
Wang and Shen's mixed foil theory, is chiefly useful for preliminary 
design. Once a few candidate hydrofoil sections have been designed, 
we must make a final selection using more exact analysis methods. 
Another point has been noted already which is peculiar to the way we have 
formulated the problem. It will be noted that because the chord line 
passes through the trailing edge of the profile in this theory, the 
cavity thickness T is a clearance between the wetted surface of the 
profile and the upper surface of the cavity. In the case of the cavity 
thickness near the nose, the distance tT is measured from the chord 
line and not the wetted surface. Therefore wT is only indirectly 
related to the clearance between the wetted surface and the upper surface 
of the cavity. As a result our study shows that occasionally one can 
obtain profiles which have cavity contours which cut the wetted surface. 


Therefore one can not guarantee from the outset that his design will be 


entirely satisfactory. As we have noted already, if one were to change 
the condition y (x, )=yT to read y (x) -n (x, )=uT and implement the design 
process accordingly, the occasional interference between cavity and 
wetted surface could be eliminated. lHowever in this case too, the 
parameter, &, would be limited in its range to certain permissible values 


which must be searched out before a design can be completed. In the 
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present somewhat simpler theory the limited range of uw makes it possible 
for one to explore the possibilities afforded by various selections of yu 


when other design parameters are held fixed. Thus he can determine from 


the resulting profile and cavity contours which of the several designs 


if any, is the most satisfactory, not only from the viewpoint of profile 
geometry but of hydrodynamic performance also. Thus we feel that the 
procedure proposed in this study offers a flexible design tool for the 
preliminary design of fully cavitating hydrofoil sections employing 

any one of the three design procedures studied. The new computer program 
for the design method has been written to operate in accordance with 
either the first, second or the third procedure. However the added 
flexibility of the third procedure which enables the designer to examine 
a variety of geometric and hydrodynamic aspects relating to his practical 
design goals by varying wu within the permissible range for systematic families 
of pressure distribution shape should make the third procedure the most 


useful of the three methods. 
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APPENDIX 


TABLES OF HYDRODYNAMIC PERFORMANCE 


NOTATION 
K cavitation number 
b cavity thickness at trailing edge 


XO location of second cavity thickness control point on 
chord line, Xo 


lift coefficient, Cc 


peak pressure location, s 
nose cavity thickness parameter, U 
lift to drag ratio 


moment coefficient, Cr 


drag coefficient, Cy 


center of pressure location along chord line, x 


angle of attack a 


cavity length parameter, Quteen™ 


point-drag solution strength, y,Q)=E/a 


fraction of lift born by prescribed pressure distribution, m 
strength of leading edge or "flat plate" solution, A 
cavity closure singularity strength, B 


value of a-Cy where C_ is the first Fourier coefficient for 
the camber function 
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z - PLANE 


Figure 2 - Mathematical Representation of the Hydrofoil and 
Cavity in the Complex z-Plane 
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Figure 3 - The Hydrofoil and its Cavity in the Vv and @-Planes 
After the Application of Conformal Mappings 
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(a) NOSE-LOADED PROFILE (b) TAIL-LOADED PROFILE 
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Figure 4 - Parametric Representations of Rectangular Pressure 
Distributions for Nose-Loaded (a) and Tail-Loaded 
(b) Hydrofoil Sections 
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Figure 5 - Limiting Values of C,/T for Rectangular Pressure 
Distributions at Zero Cavitation Number 
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Figure 6 - Permissible Values of yi at Zero Cavitation Number 
and Lift Parameter m for Rectangular Pressure 
Distributions and Three Nose Control Points, Xo 
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Figure 


8 - Semi-Elliptical Pressure Distributions for Nose and 


Tail Loadings 
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Figure 16 - Three C, Intervals which Covern the Computation 
of the Permissible Range of i 
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Figure 17 - Effect of Permissible Values of uv on L/D for 
Semi-Elliptical Pressure Distributions at K=0 


and Various C 
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Figure 18 - Effect of Permissible Values of u on L/D 
for Semi-Elliptical Pressure Distributions 


at K=0.2 and Various Cy 
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Figure 19 - Effect of Permissible Values of j) on L/D for Semi- 
Elliptical Pressure Distributions at Various 


Cavitation Numbers and Fixed Cr 
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Fixed Peak Pressure Location and Ranges of K and 


i 


5 July 1977 
BRP: JFsjep 


= 
—! 
o 
K- 
<< 
O& 
© 
= 
Qc 
7 
o) 
a 
rm 
—) 


— 
oS 


Se en a 
0.08 0.10 0.12 0.14 0.16 
LIFT COEFFICIENT, C 


Figure 21 - Effect of Permissible Values of ») on L/D for 
Various Cavity Thicknesses and Lift Coefficients 
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